Wave Propagation in the Cantor-Set Media: 
Chaos from Fractal 
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Propagation of waves in the Cantor-set media is investigated by a renormalization-group type 
method. We find fixed points for complete reflection, T* — 0, and for complete transmission, 
T* — 1. In addition, the wave numbers for which transmission coefficients show chaotic behaviors are 
reported. The results obtained are for optical waves, and they can be tested in optical experiments. 
Our method could be applied to any wave propagation through the Cantor set. 

PACS numbers: 42.25.-p, 61.44.-n, 46.65.+g, 71.55.Jv 
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Localization of the electronic states due to disorder is 
an active field in condensed matter physics. It has been 
recognized that localization occurs not only in disordered 
systems but also in quasiperiodic systems in one dimen- 
sion [l|, [3, [I[ . While the localization of states was origi- 
nally regarded as an electronic problem, it was later rec- 
ognized that the phenomenon is essentially a consequence 
of the wave nature of electronic states. Therefore, such 
localization can be expected for any wave phenomenon. 
An optical experiment with the Fibonacci multilayer was 
proposed^, and the corresponding experiments using di- 
electric multilayer stacks of SiC>2 and TiC>2 thin films 
were reported^, @. The transmission coefficients and 
scaling properties well agree with the theory and can 
be considered as an experimental evidence for localiza- 
tion of electromagnetic waves. Recently, an experiment 
for propagation of electromagnetic waves in a three- 
dimensional fractal cavity was reported]?]], where a con- 
finement of electromagnetic waves in the fractal structure 
was observed. 

Propagation of light in the Cantor media was numeri- 
cally studied by Konotop et al.[8[ and Sun and Jaggard 

a, later by Bertolotti et al. with transfer matrices [PCI 
, and recently by Yamanaka and Kohmoto [l2j and 
Hatano[lH again with transfer matrices. Comparison of 
the experimental results and numerical results was per- 
formed by Lavrinenko et al.[14j. 

We study wave propagation in the Cantor-set media by 
a renomalization-group type method[15j|. For specific val- 
ues of the wave numbers, transmission coefficients show 
chaotic behaviors governed by the logistic map. This 
exotic behavior could be observed in an optical exper- 
iment. In this system, the onc-dimentional theory is 
strictly valid. In addition, it is feasible to construct the 
system accurately, and the parameters can be precisely 
controlled and measured. 

The procedure of constructing the Cantor set begins 
with a line segment with unit length. We regard this as 
substrate A. The line segment is divided into three parts. 
The left and the right segments are substrate A, each of 
which has length 1/3. The middle part, which has length 
1/3, is removed. We regard the removed part as substrate 
B. The procedure is repeated for each of line segments A. 
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FIG. 1: The initial state j = and the first three generations 
in the construction of the Cantor set. 
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FIG. 2: Electromagnetic wave propagation through the Can- 
tor sequence Cj. 



We call j-th generation of the Cantor sequence as Cj . For 
Cj we have a set of 2 J line segments of substrate A, each 
of which has length 1/3- 7 . The Cantor set is obtained if 
this procedure is repeated infinite times. It is self-similar 
and has the fractal dimension log 2/ log 3. Constructions 
of first few generations are shown in FigJT] In the follow- 
ing, we suppose that the indices of refraction of A and 
B are ua and Ub, and the magnetic permeabilities of A 
and B are /za and ^b, respectively. 

We consider wave propagation through Cj . See FigfJ] 
For simplicity, suppose that the incident light is normal 
and the polarization is perpendicular to the plane of the 
light path. Here represents the incident light, and 



E^' and E^' represent reflected and transmitting light, 
respectively. We also consider light coming from right 

(2) 

E^'. By setting 

E+ = E<»+E<- 2 \ = -EW)/i, (1) 
the wave propagation is given by0, [HJ 
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FIG. 3: Transmission coefficient 7} as a function of the wave 
number k for C2 (top) and C3 (bottom), (jia = 2, tib = 1 and 
n — 2.) The range of the wave number k for Cj is < k < 3 J/ k. 
Note that the horizontal axes are rescaled. 



where Mj satisfies 



M j+ m=M j {\)T(^)M j {\) t (3) 



with the initial condition 

M (k) = T BA T{n A k)T AB . 



(4) 



Here k is a wave number in vacuum, and the matrices 
7ab and T BA represent light propagation across interfaces 
A <— B and B <— A respectively, 



1 
1/n 



n=^x^. (5) 
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The matrix T(5) represents light propagation within a 
layer, 



T{8) 



cos 6 — sin 5 
sin S cos <5 



(6) 



For a layer of type A and B with thicknesses d, the phase 
S is given by S = riA^rf and <5 = n B kd, respectively. 
Equation ([3]) can be regarded as a renormalization-group 
transformation. One can easily show that 



det Mi 



(7) 
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By putting ; = 0, the transmission coefficient Tj of 
Cj is given by 



wave numbers by three each time one increases a genera- 
tion. Thus we study wave propagation of wave numbers 
3 J fcofC.,. 

In Cj, all A's have the same length 1/3- 7 ', thus the 
propagation in A's are given by the unique trans- 
fer matrix T(n A k). On the other hand, B's have 
lengths 1/3 J , 1/3 J_1 , • • • and 1/3, thus in general they 
give j different transfer matrices, which are given by 
T{n B k), T(3n B k), • ■ • , and Tffi n B k). However, if k 
is given by 



run 
~ 35 ' 



(9) 



with integers m and q, T(8)'s in large B's are equal to 
T(to7t) . Now the map ([3]) for j > q is written as 



A/}+i = Nf , 

where 

Afj =T(mir)M 3 . 
From ©, Afj satisfies 



det A/} = 1. 



(10) 



(11) 



(12) 



To solve pH|) . we rewrite J\fj in terms of the Pauli matri- 
ces <Ji (i = 1.2.3), 

Afj = Xjl + ttj-a, dj = {af ] , laf , of 5 ) , (13) 

where Xj and (i — 1,2,3) are real. Comparing the 
coefficients of 1 and a in both side of flUl). we have 



Xj+l =Xj +dp 



and 



dj+i = 2xjdj. 
In addition, (fT2"]) leads to 



(14) 
(15) 



= 1. 



(16) 



From (fT4f and (fTB]) . we have the map for Xj = iTr Mj 



x j+ i = 2Xj-l, (j > q). 



(17) 



This is the logistic map which maps the interval [— 1, 1] 
onto itself and others to +00. The bounded Xj's are given 
by 



Xj = cos[2 J 9 (cos 1 Xq)}, 



(18) 



T, = 



\Mj\ 2 +2> 



(8) 



where 1 j | 2 is the sum of the squares of each matrix el- 
ement of Mj . An example of transmission coefficients as 
a function of wave numbers is shown in Fig[3]for C2 and 
C3. They show a scaling behavior when one multiplies 



and have chaotic behaviors with the Lyapunov exponent 
In 2. From (|15|) . we also find the following constants of 
motion A and rj 
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^2)= A > -lT) =r l Ci>?)- 
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FIG. 4: Transmission coefficient Tj as a function of n for 
nefc = 7r/3. Iteration number is j = 5 (top), j = 10 (center), 
and j = 100 (bottom). 



FIG. 5: Transmission coefficient Tj as a function of n for 
nsk — 7r/4. Iteration number is j = 5 (top), j = 10 (center), 
and j = 100 (bottom). 



In terms of them, the transmission coefficient ((8} is given 
by 



T, 



l-(A 2 + 77 2 ) 
1 - (A 2 + r?)x] 



(20) 



If the constants of motion A and r\ satisfy < X 2 + rj 2 < 1, 



we have 



< 1. Then, the complete transmission Tj = 
1 is achieved when Xj = ±1, and we have the minimal 
transmission 



when 



T mi „ = l-(A 2 +r? 2 ), 
0. On the other hand, if A 2 



T 



(21) 
> 1 we 



have \xj\ > 1. Then escapes to infinity and the trans- 
mission coefficient flows into the fixed point for complete 
reflection T* — 0, as seen from (|2U)) . We show the trans- 
mission coefficient Tj as a function of n for nsA; = 7r/3 
in FigdJ Since the complete reflections and the complete 
transmissions are nearby in n, the Cantor media could 
be used as a fast swiching devices. 

We note that similar chaotic behaviors also appear for 
the initial wave numbers 



1/2 



3" 



(22) 



with integers m and q. The matrix T(5) in large B's is 
now T((m+l/2)7r), and the map for Mj is Afj+i = — Nf- 
The trace of Mj follows the logistic map, again. 




FIG. 6: Transmission coefficient Tj as a function of n for 
nefc = 7r/3 + 0.0001. Iteration number is j — 5 (top), j = 10 
(center), and j — 100 (bottom). 
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For initial wave numbers other than ([9]) or (|22|) , our nu- 
merical calculations suggest that light reflects completely 
or transmits completely by the Cantor-set media. An ex- 
ample of our numerical results are shown in Fig)5] The 
transmission coefficients tend to flow into the fixed point 
T* = or T* = 1 as the generation increases. 

When initial wave numbers are near a value ^ or 
(122)1 . the transmission coefficients show chaotic behav- 
iors in first few generations. For example, the transmis- 
sion coefficients Tj for n^k — tt/3 + O.OOOf are shown 
in Fig[51 In the fifth generation (j = 5), the transmis- 
sion coefficients show chaotic behaviors indistinguishable 
from the case for nefc = 7r/3 (Fig|4]). However, in the 
tenth generation (j — 10), the transmission coefficients 



show different behaviors from the case for nefc = 7r/3, 
and in the hundredth generation (j = 100), they tend to 
flow into the fixed point T* = or T* = 1. Our numer- 
ical results again suggest that light reflects completely 
or transmits completely by the Cantor-set media. This 
behavior should be observed experimentally. 

We can also study wave propagation for other types 
of Cantor-set media by using similar recursion relations 
as (J3)). Chaotic behaviors governed by the logistic map 
occur in the Cantor-set media, if the length of each A 
becomes 1/a each time one increases a generation, and a 
is an integer. Details will be given in a separated paper 

We acknowledge useful discussions with M. Yamanaka. 
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